Stability of spin-orbit coupled Fermi gases with population imbalance 
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We use the self-consistent mean-field theory to analyze the effects of Rashba-type spin-orbit coupling (SOC) 
on the ground-state phase diagram of population-imbalanced Fermi gases throughout the BCS-BEC evolution. 
We find that the SOC and population imbalance are counteracting, and that this competition tends to stabilize 
the uniform superfluid phase against the phase separation. However, we also show that the SOC stabilizes 
(destabilizes) the uniform superfluid phase against the normal phase for low (high) population imbalances. In 
addition, we find topological quantum phase transitions associated with the appearance of momentum space 
regions with zero quasiparticle energies, and study their signatures in the momentum distribution. 

PACS numbers: 05.30.Fk, 03.75.Ss, 03.75.Hh 



Introduction. The recent realization of synthetic gauge 
fields with neutral bosonic atoms QJJ] , evidently seen from the 
appearance of vortices in a BEC, has sparked a new wave of 
theoretical interest in the cold atom community. This novel 
technique uses a spatially-dependent optical coupling between 
internal states of the atoms, and can be used to engineer more 
complicated gauge fields by dressing two atomic spin states 
with a pair of lasers. For instance, it has recently been used 
to create and study the effects of SOC in a neutral atomic 
BEC with equal Rashba and Dresselhaus strengths [2]. Since 
this method is equally applicable for neutral fermionic atoms, 
and given that the coupling between a quantum particle's spin 
and its momentum is crucial for the topological insulators and 
quantum spin Hall states, which have recently received im- 
mense interest in the condensed matter community fl, it may 
allow for the realization of topologically nontrivial states in 
atomic systems with possibly a broad interest in the physics 
community 

Motivated by the recent success in realization of the SOC 
Bose gases [2], and by a practical proposal for generating a 
SOC in 40 K atoms (70, effects of the SOC have recently been 
studied for the two-component Fermi gases: (i) the two-body 
problem exactly [8], and (ii) the many-body problem in the 
BCS mean-field approximation raJfjllS Wl ill . It has been found 
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that the increased density of states due to the SOC plays a cru- 
cial role for both problems. In particular, for the two-body 
problem, this gives rise to a two-body bound state even on the 
BCS side (a s < 0) of a resonance |8J]. For the many -body 
problem, the increased density of states favors the pairing 
so significantly that increasing the SOC, while the scattering 
length is held fixed, eventually induces a BCS-BEC crossover 
even for a weakly-interacting system when a s — > CP [9- 
1 ill . In addition, the SOC leads to an anisotropic superfluid, 
the signatures of which could be observed in the momentum 
distribution or the single-particle spectral function for suffi- 
ciently strong SOCs lUHl . 

In this paper, we study the competition between the SOC 
and population imbalance on the ground-state phase diagram 
of two-component Fermi gases across a Feshbach resonance, 
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FIG. 1 : (Color online) The ground-state phase diagrams are shown 
as a function of population imbalance P = (JVf — N{)/N and scat- 
tering parameter l/(kF<i s ), where the SOC parameter ma/ltF is set 
to 0.005 in (a), 0.05 in (b), 0.15 in (c), and 0.25 in (d). We show 
normal (N), phase separation (PS), gapped superfluid (SF) and gap- 
less superfluid (gSF) phases. The dashed blue and dotted green lines 
separate topologically distinct gSF regions, and the trivial SF phase 
resides at and around the P = line. 



i.e. throughout the BCS-BEC evolution. Our main results 
are shown in Fig. Q] and they are as follows. In the absence 
of a SOC, the phase diagram of population-imbalanced mix- 
tures is well-studied in the literature IU2I1 . and it mainly in- 
volves normal (N), phase separation (PS), and topologically 
distinct gapless superfluid (gSF) as well as the trivial gapped 
superfluid (SF) phases. In particular, at and around unitarity 
(\a s \ — > oo), the system changes from SF to PS and then to 
N as a function of increasing population imbalance. In this 
paper, we show that the SOC and population imbalance are 
counteracting, and that this competition always tends to stabi- 
lize the gSF phase against the PS. However, we also show that 
the SOC stabilizes (destabilizes) the gSF phase against the N 
phase for low (high) population imbalances. In addition, since 
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the SOC stabilizes the gSF phase for a very large parameter 
region at and around the unitarity, where most experiments are 
conducted [ 14], it may allow for a possible realization of the 
gSF phase for the first time with cold atoms. 

Mean-field Hamiltonian. To obtain these results, we use the 
mean-field Hamiltonian (in units of ft = 1 = fc#) 
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which is defined up to the constant C = (1/2) £ kj(T + 

|A| 2 /g. Here, V k = [a kt , a]^, a_ k , t , a_ k ,j.] denotes the 

fermionic operators collectively, where a k (a k , CT ) creates 
(annihilates) a spin-cr fermion with momentum k, £ kcr = 
£- k .<r = e k (T — [i a has the inversion symmetry with e k cr = 
k 2 /(2m a ) the kinetic energy, /i^ the chemical potential, and 

k = ^Jk 2 + k 2 + k 2 . Here, A = <?(a k ^a_ k ^) i s tne mean- 
field order parameter, where g > is the strength of the attrac- 
tive particle-particle interaction which is assumed to be local, 
and (• • • ) is a thermal average. While we keep the formalism 
quite general including the possibility of mass-imbalanced 
fermion mixtures, we present our numerical results only for 
mass-balanced mixtures. In addition, we consider only a 
Rashba-type SOC, i.e. S k = — S-k — a(k y — ik x ), where 
a > is its strength, but the extension of our work to 
anisotropic or more complicated SOCs is straightforward. 

Self-consistency equations. Transforming the mean-field 
Hamiltonian to the helicity basis, and following the usual 
procedure, i.e. 9f2/<9|A| = for the order parameter and 
N+ + sNi = —dil/dfig for the number equations where 
f2 is the mean-field thermodynamic potential, s = ± and 
Ms = (mt + S AH) /2> we obtain the self-consistency equations 
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Here, T is the temperature and E 2 . s = C k + + C k 



S k 



2syl k gives the quasiparticle excitation spec- 



tram Si, where A k - ^,_(^,+ + |A| 2 ) + |S k | 2 £ 2 + , 

Ck. s = £k,s - Ms with e k , s = (ek,f + se k ,j.)/2 = k 2 /(2m s ) 
and m s = 2m^mi/(mi + sm-j-). Note that m + is twice 
the reduced mass of j" and J, particles, and m_ — > oo for 
mass-balanced (m-f- = mj.) mixtures. In Eqs. d2} - ©, 
the derivatives of the quasiparticle energies are given by 
&E k , s /9|A| = (1 + sCk,-Mt) I A l for the order param- 
eter, dE ktS /dn+ = -|l + + |S k | 2 )M k ]£ k> _/£ k , s 
for the average chemical potential and dE k:S /dfi- = — [1 + 



s(£ k ,+ + |A| 2 )/A k ]£ k ,-/E k , s for the half of the chemical po- 
tential difference. 

Equations (f2]l - (0]l are the generalization of the mean- 
field expressions to the case of population- and/or mass- 
imbalanced mixtures, and they recover the known expressions 
(a) when £ k ._ = for which i? k s simplifies to E^ s = 
(£k,+ + s\S k \) 2 + |A| 2 , and (b) when \S k \ = for which 
E Ks simplifies to E 2 s = [s&,_ + Jtf~ + |A| 2 ] 2 . We 
eliminate the theoretical parameter g in favor of the exper- 
imentally relevant s-wave scattering length a s via the re- 
lation, 1/g = -m + Vl{Aira s ) + S k l/(2e k)+ ), where V 
is the volume, g can also be eliminated in favor of the 
two-body binding energy < in vacuum via the re- 
lation [8] 1/g = (1/2) £ ki , l/(2e k ,» + eth - e&), where 

\S k \ 2 is the single-particle nonin- 
teracting dispersion in the helicity basis, and eth is the thresh- 
old for the two-body scattering. For the Rashba-type SOC, 
and assuming m-f < m,[., we obtain e t h — 2a 2 m_(m_ — 

^Jrn 2 _ — m 2 + )/m + , which gives e t h = m + a 2 when m_ ^> 

m + , i.e. m-f < mj,. 

Thermodynamic stability. To construct the phase diagram, 
we solve the self-consistency equations and check the sta- 
bility of these solutions for the uniform superfluid phase us- 
ing the compressibility (or the curvature) criterion Tdft . 
This says that the compressibility matrix k(T) with elements 
K a . a '(T) = — d 2 il/ (d^a-dfia') needs to be positive definite, 
and it is directly related to the condition that the curvature of 
il with respect to | A|, i.e. 
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needs to be positive. When at least one of the eigenvalues of 
k(T), or the curvature <9 2 f2/<9|A| 2 is negative, the uniform 
mean-field solution does not correspond to a minimum of ft, 
and a nonuniform superfluid phase, e.g. a phase separation, is 
favored ifBl [Till. 

Rashba-type SOC for mass- and population-balanced mix- 
tures. To get more insight into the effects of the SOC 
in the BCS-BEC evolution, let us first consider mass- and 
population-balanced {m a = m and N a = N/2) mixtures, 
where N is the total number of fermions. This case is an- 
alytically more tractable, since the single-particle energy in 
the helicity basis simplifies to [151] £k,s = k 2 /(2m) + sak±, 

where k± = ^Jk 2 + k 2 . For instance, the bound-state equa- 
tion can be solved to obtain l/a s = \Jm 2 ct 2 — met, + 
ma\n\^—meb/{^/m 2 a 2 — meb + ma)]. In the weak SOC 
limit, when ma 2 <C eb, this expression gives eb ~ 



■i/W 



up to the leading order in a, which recov- 



ers the usual result in the a — > limit. However, in the strong 
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SOC limit, when ma 2 ^> ej,, the general expression gives 
€b ~ -(4mct 2 /e 2 )e" 2 ^ maa> \ which implies that a bound 
state exists even for a s < 0, although its energy is exponen- 
tially small [8]. This is a result of increased density of states 
D s(e) = Ek'K 6 - £ k, s ) due to the SOC HELl, where S(x) 
is the Dirac-delta function. 

In the case of noninteracting (g = or a s —> Q~) mass- and 
population-balanced Fermi gases at zero temperature (T = 0), 
where /i a = [i, calculating the total number TV = N + + N- — 
kpV/(3TT 2 ) of fermions, where N s — X^k^A* ~ £ k.s), we 
obtain /i « ep - 3ma 2 /2 up to the leading order in a when 
ma 2 <C e_F, and /i = — ma 2 / 2 + 2kp / (Zirm 2 a) when /i < 0. 
Note that we conveniently choose the energy (length) scale 
as the Fermi energy ep (momentum Uf) of the N a = N/2 
fermions. 

It has been shown that increasing the SOC for a noninteract- 
ing Fermi gas leads to a change in the Fermi surface topology, 
when the number of fermions in the +-helicity band (iV + ) 
vanishes H. This occurs when /i goes below the bottom of 
the energy band, i.e. /i = 0, or when a increases beyond 
a = [4/ (Sn^^kF /m ~ Q.Tbkp/m. In some ways, this is 
similar to the usual BCS-BEC crossover problem, where the 
quasiparticle excitation spectrum changes behavior as a func- 
tion of increasing the scattering parameter 1 / (kpa s ) at /i = 0, 
i.e. its minimum is located at a finite (zero) momenta when 
/i > (/i < 0). However, the topological transition discussed 
here is driven by increasing the SOC parameter a, and the 
origin of it can also be traced back to the change in the quasi- 
particle excitation spectrum when g is finite. For instance, for 
mass- and population-balanced mixtures, the excitation spec- 
trum simplifies to itMiKil #k, s = V( £ k,,-^) 2 + |A| 2 , 
and it also has a change of behavior at /i — 0. Having set 
up the formalism, we are now ready to discuss the competi- 
tion between normal fluidity, uniform superfluidity, and phase 
separation across a Feshbach resonance. 

Ground-state phase diagrams. In the absence of a SOC and 
at low T, it is well-established that the self-consistent solu- 
tions of Eqs. (O - dU are sufficient to describe the physics of 
fermion mixtures both in the BCS and the BEC limits, and that 
these equations also capture qualitatively the correct physics 
in the entire BCS-BEC evolution JT^Hj]. Hoping that the 
mean-field formalism remains sufficient in the presence of a 
SOC, here we analyze only the ground-state phase diagram of 
population-imbalancedbut mass-balanced mixtures as a func- 
tion of both the SOC and scattering parameters. 

There are typically three phases in our phase diagrams. 
While the normal (N) phase is characterized by A = 0, the 
uniform superfluid (SF) and nonuniform superfluid, e.g. phase 
separation (PS), are characterized by <9 2 il/<9|A| 2 > and 
d 2 il/d\A\ 2 < 0, respectively, when A ^ 0. Furthermore, in 
addition to the topologically trivial gapped SF phase, the gap- 
less SF (gSF) phase can also be distinguished by two topo- 
logically distinct regions, depending on the momentum-space 
topology of their quasiparticle excitation spectrum (see be- 
low). 

In Fig. Q] the phase diagrams are shown as a function of 
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FIG. 2: (Color online) The ground-state phase diagrams are shown 
as a function of P = (TV-f — N±)/N and a, where l/(fcfa s ) is set 
to —0.5 in (a), in (b), 0.5 in (c), and 1.5 in (d). The labels are 
described both in Fig.fTJand in the text. 



population imbalance P = (Nf — N±)/N and scattering pa- 
rameter l/(kpa s ) for four different a values. Comparing 
these results with the a — > limit lfl3ll . it is clearly seen that 
the SOC and population imbalance are counteracting. On one 
hand, this competition always tends to stabilize the gSF phase 
against the PS, and therefore, at any given P, the system even- 
tually becomes a gSF by increasing a, no matter how small 
l/(kpa s ) is. This is best seen in Fig. [2] where the phase di- 
agrams are shown as a function of P and a for four different 
l/(kFa s ) values. On the other hand, we find that while the 
SOC stabilizes the gSF phase against the N phase for low P 
due to increased density of states lfl7ll . it destabilizes the gSF 
phase against the N phase for high P. 

In both figures, the dashed blue and dotted green lines are 
obtained from the conditions |A| 2 = —/.i-ffii and |A| = 
(see below), respectively, and they separate the topologically 
distinct gSF regions. The trivial gapped SF phase resides at 
and around the P — lines. As can be inferred from Fig. |2] 
the dashed blue line makes a dip as l/(k Fd s ) S 0, the tip of 
which eventually touches the P = line at a w 0.75fcp /m, 
consistent with our analysis above for the topological transi- 
tion of a noninteracting (a s — > CP) system. In Figs. Q] and [2] 
we also show that the SOC stabilizes the gSF for a very large 
parameter region, at and around the unitarity, which is nor- 
mally unstable against PS when a = 0, allowing for a possi- 
ble realization with cold atoms for the first time, as discussed 
next in great details. 

Topological phase transition. Population imbalance is 
achieved when either £'k.+ (for iVj, > Nf) or -Ek,- (for 
N^ > N±) has zeros in some places of k-space. Let us as- 
sume Nf > Ni without loosing generality, for which i?k,+ is 
always gapped. Depending on the number of zeros of -Ek,- 
(zero energy surfaces in k-space), there are two topologically 
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distinct gSF phases: gSF(I) where i?k,- has two, and gSF(II) 
where E^^ has four zeros. The zeros of -Ek.- can be found 
by imposing the condition E^ + E"^_ = (£k,t£k,4 + 1^1 2 — 
|Sk| 2 ) 2 + 4|A| 2 |S' k | 2 = 0, indicating that both |5 k | = and 
Ck.tCk,^ + | A| 2 = needs to be satisfied. 

For the Rashba-type SOC that we consider in this pa- 
per, the zeros occur when k± = and at real k z mo- 
menta, k\ s = B + + s \J B 2 ^ — 4m-|-TOj_| A| 2 , provided that 

|A| 2 < | J B_| 2 /(4m t m i ) for B+ > 0, and |A| 2 < -fiffj,± 
for B + < 0. Here, B s — m^/if + sm^fi^. Note that the con- 
ditions on k z coincide with those obtained for k in the case 
of population-imbalanced mixtures without the SOC where 
i?k.s is isotropic. The topologically trivial SF phase corre- 
sponds to the case where both £"k.+ and -Ek,- have no zeros 
and are always gapped. The transition from gSF(II) to gSF(I) 
occurs when \k Xt -\ —> 0, indicating a change in topology in 
the lowest quasiparticle band, similar to the Lifshitz transition 
in ordinary metals and nodal (non-s-wave) superfiuids. 




k Zt+ = 1.00k F in Fig.Oa), and k z _ = and k Z:+ = 0.97fc F 
in Fig. Ob), in perfect agreement with our analysis above. 
We also see that a major redistribution occurs for the minor- 
ity species (nk,j,) at the gSF(II) to gSF(I) transition bound- 
ary, where the sharp peak that is present near the origin van- 
ishes abrubtly. Although this topological transition is quan- 
tum (T = 0) in nature, signatures of it should still be observed 
at finite T within the quantum critical region, where the nk, ff 
are smeared out due to thermal effects. Although the primary 
signature of this topological transition is seen in nt^, single- 
particle spectral function 111 811 as well as some thermodynamic 
quantities such as the atomic compressibility would also show 
an anomaly at the transition boundary. 

Conclusions. In summary, we analyzed the effects of SOC 
on the ground-state phase diagram of population-imbalanced 
Fermi gases throughout the BCS-BEC evolution. We found 
that the SOC and population imbalance are counteracting, and 
that this competition always tends to stabilize the gSF phase 
against the PS. In contrast, while the SOC stabilizes the gSF 
phase against the N phase for low population imbalances, it 
destabilizes the gSF phase against the N phase for high popu- 
lation imbalances. In addition, we found topological quantum 
phase transitions associated with the appearance of momen- 
tum space regions with zero quasiparticle energies, and stud- 
ied their signatures in the momentum distribution. We hope 
that our work will motivate further research in this direction, 
since the SOC stabilizes the gSF phase for a very large param- 
eter region at and around the unitarity, allowing for a possible 
realization of the gSF phase for the first time with cold atoms. 
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FIG. 3: (Color online) Typical momentum distributions nk, CT are 
shown as a function of k± — y/k% + fc| and k z , where we set 
P = 0.5 and l/(kpa 3 ) = 0, and vary the SOC parameter: (a) 
a = 0.275fc F /m (for which |A| = 0.463e F , = 0.550eF 
and fi- = 0.653£f), and (b) a = 0.350kp /m (for which |A| = 
0.465eF, fi+ = 0.427ef and /Lt- = 0.688£f), corresponding to 
gSF(II) and gSF(I) phases, respectively. 

The topological transition here is unique, because it in- 
volves an s-wave superfluid, and could be potentially ob- 
served for the first time through the measurement of the mo- 
mentum distributions n^. a of \ and I fermions, both of which 
are readily available from Eqs. (|3} and (@). For instance, we 
illustrate the typical T = distribution of a gSF(II) phase 
in Fig. Oa), and of a gSF(I) phase in Fig. [3jb). In these 
figures, we note that n^^ is anisotropic in k space, which 
follows from the anisotropic structure of E^, s . For k-space 
regions where k± = and k z ,- < \k z \ < k z , + , the corre- 
sponding distributions are exactly nk,t = 1 an d n k.l — 0. 
Here, k Z:S are approximately found to be = 0.30fcF and 
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